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Summary — Finite-time thermodynamics optimisation of an endo- exo-irreversible Stirling Motor. This work deals 
with the finite time thermodynamic optimization of a Stirling engine. An endo- and exo-irreversible cycle is considered. 
The internal irreversibilities are due to the internal conductance of the plant and to the non-adiabatic regenerator. Also, 
the general irreversibilities of the non-quasistatic cycle are taken into account. The external irreversibilities are due to the 
linear interactions between the heat sources and the working fluids at finite temperature gaps. 

A sensitivity analysis was performed in order to optimize the operation conditions leading to maximum power output. 
The results obtained are in good agreement with the experimental data. 


1 ■ INTRODUCTION 


Bien que l’interet pour l’ensemble des problemes 
du moteur de Stirling ait connu des fluctua¬ 
tions, tenant compte des plus recentes realisations 
[1-4], on remarque neanmoins une preoccupation 
soutenue pour leur perfectionnement. 

Tout en considerant les difficultes liees a ce per¬ 
fectionnement des machines de Stirling comme une 
motivation de plus en vue du developpement, on 
considere aussi qu’il faut accorder une attention 
particuliere a la recherche theorique et experimen- 
tale dans le domaine d’application de ces machines 
thermiques constituant une alternative pour resou- 
dre les problemes energetiques et ecologiques de la 
society. 

En suivant l’optimisation des systemes 
thermoenergetiques nucleaires, dans le cas par- 
ticulier endoreversible, Novikov [5] a obtenu 
en 1957, probablement le premier, l’expression 
connue aujourd’hui sous le nom de The Nice 

Radical jV = 1 - y/T 0 /T < r) m ax = 1 - (T 0 /T)j , qui 

represente la valeur du rendement thermique as- 
socie au regime de production d’une puissance 
mecanique maximale. Plus tard, ce resultat a ete 


* Correspondance et tires a part 


confirme, notamment grace aux travaux de Curzon 
et Ahlbom [6]. Maintenant, ce resultat est bien 
accepte dans la thermodynamique en temps fini. 

Le trait specifique des modeles d’analyse bases 
sur la thermodynamique en temps fini consiste dans 
la prise en consideration et dans l’optimisation des 
differences finies de certains parametres intensifs 
(temps A t , temperature AT, pression Ap, etc). 
Ceux-ci sont lies aux dissipations internes et 
externes qui apparaissent dans le fonctionnement 
des systemes considers, optimisation ayant pour 
but la maximisation des performances. 

Ce travail est un prolongement des modeles 
d’analyse du moteur de Stirling a cycle ideal 
endo-reversible mais exo-irreversible developpes sur 
la base de la thermodynamique en temps fini 
[7-11]. II propose une approche plus complete de 
l’optimisation des conditions de fonctionnement du 
moteur pour produire une puissance mecanique 
maximale. 


2 a MODELISATION DU CYCLE POUR 
LE MOTEUR DE STIRLING 


Les hypotheses qui constituent les fondements 
du modele d’analyse theorique du moteur de 
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Stirling endo- et exo-irreversible considere sont les 
suivantes. 


a) On considere le cycle ideal du moteur de Stir¬ 
ling ( fig 1), forme de deux processus isothermiques 
( T c = Ti, 2 = cte; Tp = T 3 , 4 = cte) et de deux pro¬ 
cessus isochores (V m — V 2 ,z = cte ; Vm = Vi l4 = cte), 
avec les parametres adimensionnels suivants : 

— e = rapport volumetrique de compression 

» m 

isothermique; 

— f3 = rapport des pressions extremes du 

Pm 

cycle; 


Ti_ 

T r 


rapport des temperatures moyennes 


des sources de chaleur. 



b) Par rapport au cycle ideal de Stirling endo- 
et exo-reversible, dans ce modele, on introduit 
les sources d’irreversibilite interne et exteme 
suivantes : des pertes de chaleur et des sources 
d’irreversibilite interne. 

— Des pertes Q p dues a la conductance thermique 
interne entre la partie chaude, /, et la partie 
froide, R, du moteur. Ces pertes sont directement 
proportionnelles a la difference de temperature 
entre les deux extremites de la machine [ 12 ], et 
le flux de chaleur perdu aura done l’expression 
suivante : 

Q p = K p (Td — Tc) (21) 


l’agent cycle a son passage vers l’espace de com¬ 
pression Q 4 ,i et la chaleur saisie du regenerateur 
par 1 ’agent cycle a son passage vers l’espace de 
detente Q 2 , 2 >- Dans l’hypothese du gaz ideal, le 
coefficient adiabatique du reg 6 nerateur est donne 
l’expression: 


Q 2 ,2' T 2 , — Tc 

_ _ T D -Tc 


( 2 . 2 ) 


Pour le cas ideal de l’isolation parfaite du rege¬ 
nerateur, H a — 1 . 

— Des pertes dues a l’irreversibilite du cycle 
thermodynamique determinees par les processus 
non quasi statiques du cycle ; ainsi, de liquation 


de Clausius pour les cycles irreversibles 
il resulte : 


/* 


< 0, 


(Qd+QvA (\Qc\ IQp g l\ 

V To T 2 , 3 ) \T C + T AA , J 


< 0 


(2.3) 


ou Q 2 /,3 represente la puissance thermique 
consomm 6 e en plus, pour compenser le flux de 
chaleur Q p 9 = Q 44 < perdu au regenerateur non- 
adiabatique. 

Dans la relation (2.3), les temperatures moyen¬ 
nes des processus (2'-3) et (4-4'), dans l’hypothese 
du gaz ideal, sont exprimees par la relation : 


^ 2',3 — 


Q2’ ,3 

A$2',3 


T 4 , 4 ' 


Q4.4' 

AS 4>4 / 


Tp — T 2 , 
\n(T D /T 2 ,) 


(2.4) 


Pour prendre en consideration rirreversibilite 
du cycle, on emploie un parametre global <p. Ce 
parametre est determine [13] par l’irreversibilite de 
l’echange de chaleur avec les sources : 

V? = ASr/ASi ^ 0 (2.5) 

ou A Sr et AS 1 sont respectivement les echanges 
d’entropie du cote froid et du cote chaud de la 
machine. Dans le cas du cycle endoreversible ip = 1. 
Ainsi, a l’aide de ce parametre, le bilan entropique 
(relation 2.3) devient: 


\Qc\ IQp g l\ 
Tc T aa , ) 


- V 


Qp 

Tp 


+ 


Q2',3 \ 
T 2 ',3 J 


— 0 


(2.6) 


c) Des sources d’irreversibilite exteme. 

— L’echange de chaleur agent cycle—source 
chaude (le rechauffeur du moteur etant a la 
temperature moyenne T}) a lieu a difference finie 
de temperature (AT/ = T/ - Tp). Ainsi, la puis¬ 
sance thermique obtenue dans le rechauffeur aura 
l’expression : 


ou K p (W.K -1 ) > 0, represente la conductance 
thermique de perte ; elle est fortement influencee 
par les solutions constractives adoptees, a la limite, 
dans le cas ideal de l’isolement parfait K p = 0. 

- Des pertes de regenerateur, R g , non- 
adiabatique par rapport au milieu exterieur; ces 
pertes sont pretisees par le coefficient adiabatique 
du regenerateur, /x a < 1 , defini comme le rapport 
entre la chaleur cedee dans le regenerateur par 


Qi = kjAi{Ti — Tp) (2.7) 

- L’echange de chaleur agent cycle-source froide 
(le refrigerant du moteur etant a la temperature 
moyenne Tr) a lieu aussi a difference finie de 
temperature (A Tr = T c - Tr). Et done la puis¬ 
sance thermique cedee dans le refrigerateur aura 
l’expression : 

Qr = IcrArITc — Tr) (2-8) 
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Dans les relations (2.7) et (2.8), on a considere 
le cas des interactions thermiques lineaires [7], Ici, 
h et Ai represented respectivement le coefficient 
global et la surface efficace de l’echangeur de 
chaleur i considere. Leur produit represente la 
conductance thermique de l’echangeur de chaleur i 
considere (Ki = kiAi). 


3 b indices de performance 

DU MOTEUR DE STIRLING 


Les relations (2.1)—(2.8) representent la base du 
modele propose pour l’analyse du moteur de Stirling 
endo- et exo-irreversible. Pour faciliter l’interpreta- 
tion et la comparaison des nouveaux resultats avec 
ceux publies [7-11], on va utiliser des grandeurs 
adimensionnelles. 

En considerant la necessite de compenser les 
pertes de chaleur du regenerateur non adiabatique, 
la puissance thermique re?ue par l’agent cycle a la 
source chaude sera : 

Qsc — Qd + 4 ', 3 (3.1) 

ou le flux thermique Qd, pris par l’agent cycle 
dans le processus de detente isothermique 3-4, se 
determine avec la relation : 

Qd = vuRTd ln(e) (3.2) 

et le flux thermique Q 2 , i3 , pris par l’agent cycle 
dans le processus de chauffage isochore 2'-3, 
necessaire pour compenser les pertes thermiques 
du regenerateur, a l’expression : 

Qv, 3 = rh^-j(T D -Tv) (3.3) 

Dans les relations (3.2) et (3.3), m, R et k sont 
respectivement le debit massique, la constante des 
gaz parfaits et l’exposant adiabatique de l’agent 
cycle. 

De la relation (2.2) de definition du coefficient 
adiabatique du regenerateur, on peut determiner 
l’expression de la temperature T 2 , : 

T 2 > = fiaTo + (1 — t^a)Tc (3.4) 


Tout en considerant aussi les pertes dues a 
la conductance thermique interne du moteur, le 
bilan des puissances thermiques echangees dans 
le rechauffeur permet d’exprimer la puissance 
thermique Qsc re?ue par l’agent cycle a la source 
chaude : 

Qsc = Qi — |Q P | (3.5) 

ou le flux Qi re$u dans le rechauffeur, sous forme 
adimensionnelle, est: 


Qi 


Qi 

KtTt 


AT/ 

T/ 


= x 


(3.6) 


et le flux thermique \Qp\ perdu a cause de la 
conductance thermique interne du moteur, sous 
forme adimensionnelle, a l’expression : 


1 7=) | _ \Qp\ _ Kp t d ~ Tc _ -j-p Tp - T c 
IVpl KjTj Ki Tj p Ti 


(3.7) 


Dans les relations (3.6) et (3.7), on a utilise les 
notations suivantes : 


x = 



Kp 

Kj 


(3.8) 


La puissance thermique |Qsf| cedee par l’agent 
cycle au puits froid, dans le cas ou les pertes de 
chaleur de regenerateur ont ete considerees, se 
determine avec la relation : 


IQsH = \Qc\ + |Qp S | (3.9) 


ou le flux |Qc| cede par l’agent cycle dans le 
processus isothermique de compression 1-2 est 
calcule a partir de la relation : 

\Qc\ = mRTc ln(e) (3.10) 

et le flux thermique Q* 9 cede au regenerateur par 
l’agent cycle, dans le processus de refrigeration 
isochore 4-4', et perdu par celui-ci vers le milieu 
exterieur, dans l’hypothese du gaz ideal, se calcule 
avec la relation (3.3). 

En considerant les pertes thermiques dues a la 
conductance interne du moteur, le bilan des puis¬ 
sances thermiques echangees permet d’exprimer le 
flux cede par l’agent cycle a la source froide : 

IQsf| = |Qh| - \Q P \ (3.11) 

ou le flux |Qr| cede sous forme adimensionnalis6e, 
a l’expression : 


in i jggj _ K r T r &T r __ K 
KiTi Ki T/ T r t V 


(3.12) 


Dans la relation (3.12), on a utilise les notations : 


V = 


A Tr 
Tr 



(3.13) 


Sur la base des differences relatives de tempera¬ 
ture x et y (relations 3.8 et 3.13), les temperatures 
des processus isothermiques du cycle thermodyna- 
mique endo- et exo-irreversible considere auront les 
expressions : 


T d = Ti-AT d = T/(l-x); T c = T R +AT R = T R (l+y) 

(3-14) 

En utilisant les expressions (3.14), dans la rela¬ 
tion (3.7), on obtient, sous forme adimensionnelle, 
les pertes thermiques dues a la conductance interne 
du moteur : 


\Q v \ = K p [\-x-(1 + v)/t\ (3.15) 
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Aussi, sur la base des relations (3.14), pour les 
temperatures TV (relation 3.4), TV, 3 et T 4 4 / (relation 
2.4), on obtient les expressions : 

T v = [Tfta( 1 - *) + (1 - Ma)(l + y)}T R (3.16) 

T =T = (1 ~ Ma)[r(l - X) - (1 + y)]T R 

2 ' 3 4 ’ 4 ' ln[r(l-x)/(r/i a (l-a:) + (l-p a )(l + 3 /))] 

(3.17) 

De la relation (3.3), en utilisant les relations 
(3.9), (3.11) et (3.16), pour le flux thermique Q 2 \3 il 
resulte l’expression : 

Qv .3 =- rr\y a \ -w *~ 

ln(£)fc “7( T IT 

(3.18) 

En utilisant les relations (3.5) et (3.11), le bilan 
des puissances du cycle moteur permet d’exprimer 
la puissance mecanique P foumie par le moteur de 
Stirling comme : 

P — Qsc - \Qsf\ = Qi - Qr (3.19) 


Sur la base des relations (3.6) et (3.12), la puis¬ 
sance mecanique foumie par le moteur devient 
l’expression adimensionnelle suivante : 

r= l£r, = *- ! h < 3a » 


Entre les deux variables x et y de la puissance 
mecanique, il y a la relation de liaison du bilan 
entropique pour le cycle irreversible considere 
(relation 2.6) qui, apres traitement, devient: 


(r+, = 

+ < 321 > 

ou l’on a utilise les notations : 


A(x,y) = B(x, y) 



<p , v - 1 A 

l-x C(x,y)(l - x)) 
(3.22.1) 


B(x,y) 


1 y>a 

ln ( £ )' I - 1 /("IT^- 1 ) +1 -^ 

(3.22.2) 


(l-*,)(l-ii±J) 

f( x _ _ V T l — X / _ 

ln[r(l - x)\/[rn a { 1 - x) + (1 - + j/)] 

(3.22.3) 


En divisant la relation (3.21) par le produit K 1 T 1 
et en utilisant les relations (3.6), (3.12) et (3.17), on 
obtient: 


( 


1 +y 


C(x,y) 



(pX 
1 — X 


-\-Kp 



<P 

1 — x 


C(x f y) 


1 + 1 / 


) (3-23) 


Pour le rendement thermique du moteur de 
Stirling endo- et exo-irreversible considere, sur la 
base des relations (3.5), (3.6), (3.15) et (3.20), il 
resulte l’expression suivante : 


P _ P 
Qsc Qj - |Q p | 


K 

x - y 

x-K p (l-x- 7^) 

(3.24) 


RESULTATS ET COMMENTAIRES 


Sur les figures 2 a 5, sont presentes les resultats 
obtenus pour une etude de sensibilite, basee sur 
la thermodynamique en temps fini [7], qui met en 
evidence l’influence des irreversibilites considers 
en plus du cas endoreversible mais exo-irreversible. 

Les indices de performance analyses sont la 
puissance mecanique foumie P, relation (3.20), et le 
rendement thermique du moteur P, relation (3.24). 
Ils sont fonction des parametres specifiques de perte 
et_d’irreversibilite interne (K v , <p, et y a ) et exteme 
(K, x et y). Les variables x et y, sont correiees par 
la relation (3.23), qui peut uniquement etre resolue 
par des methodes numeriques. Dans ces conditions, 
l’etude de sensibilite va aussi suivre l’influence de 
la valeur maximale de la puissance mecanique 
P max foumie par le moteur, sur les conditions 
de fonctionnement correspondantes (les valeurs 
optimales x op t et y op t des differences relatives de 
temperature), et aussi sur le rendement thermique 
optimal 77 opt, qui correspond au regime de puissance 
maximale. 

Les resultats de cette application numerique 
sont obtenus en considerant les valeurs specifiques 
suivantes pour les^moteurs de Stirling : r = 4, 
£ = 1,5, k= 1,4 et K = 100. 

L’influence de l’irreversibilite du cycle non quasi 
statique, consideree par l’intermediaire du pa- 
rametre global , sur la puissance P et sur le ren¬ 
dement 77 est presentee sur la figure 2. Ainsi, pour 
1 p passant de 1 (cas ideal) a 2 , on constate une 
diminution des performances (Pmax baisse de 60 % 
et rjopt baisse de 42 %) dans les conditions ou x opt 
baisse de 0,50 a 0,28 et y opt = 0,01 = cte. 

La figure 3 pr6sente l’influence des pertes 
de chaleur dues a la conductance interne de 
la machine, exprimees par le parametre A p . A 
l’augmentation de K v de 0 (cas ideal) a 2 correspond 
une diminution de la puissance mecanique P 
foumie, dans les conditions ou le rendement 77 
est constant. En ce qui conceme l’influence sur 
le regime de puissance maximale, on constate que 
Pmax et rjopt baissent respectivement de 48 % et de 
42 %, dans les conditions ou x opt et y op t passent 
respectivement de 0,50 a 0,64 et de 0,01 a 0,02. 
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L’effet des pertes de chaleur dues au regenera- 
teur non adiabatique, introduites avec le parametre 
p a , sur les indices de performance du moteur de 
Stirling, est presente sur la figure 4. Ainsi, on peut 
remarquer que lors de la diminution de p a de 1 
(cas ideal) a 0,80^ les indices de performance bais- 
sent egalement ( P ma x baisse de 37 % et 7/ opt baisse 
de 44 %). Sur la figure 5, sont presentees les in¬ 
fluences des pertes de chaleur et des irreversibilit6s 
internes sur les indices de performance du moteur 
P et tj. Ainsi, a partir des courbes (a) pourje cas 
endo-reversible mais exo-irreversible (p = 1, K p = 0 
et /in = 1), sont graduellement introduites, par su¬ 
perposition, les sources d’irreversibilites internes 
considerees. D’abord, les courbes (b) prennent en 
compte l’irr6versibilite mteme due au cycle non 
quasi statique (<p = 1,1, K p = 0 et /i a = 1), alors que 
les courbes (c) consi derent en plus les pertes dues 



Fig 2. Influence du parametre <p sur les performances P 

et r/. 



Fig 3. Influence du parametre K v sur les performances P 

et t). 


a_la conductance interne de la machine ip - 1,1, 
K p = 0,1 et p a = 1)- Et enfin, les courbes (d) tien- 
nent aussi compte des pertes dues au regenerateur 
non adiabatique ip = 1,1, K v — 0,1 et p a — 0,_9). 
Pour les valeurs considerees des parametres p>, K v 
et /i 0 , on constate ime influence diff6rente pour le 
regime de puissance maximale. Ainsi, dans le cas 
(b), il resulte une diminution respective pour P max 
et Tjopt de 9,3 % et de 5 % ; dans le cas (c), il resulte 
une diminution respective de 5,8 % et 4,6 %, et dans 
le cas (d), il resulte respectivement une diminution 
de 29 % et de 32 %. On remarque que la perte due 
au regenerateur non adiabatique a l’influence la 
plus forte. Ainsi, en comparaison avec le cas en- 
doreversible [7], dans le cas de la superposition des 
irreversibility considerees en plus, on constate une 
diminution des performances d’environ 40 % pour 
le regime de puissance maximale P ma x et j? opt . 



Fig 4. Influence du parametre p a sur les performances P 

et rj. 



Fig 5. Influence des parametres p, K v et p a sur les 
performances P et rj. 
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5 . CONCLUSION 


Ce travail represente une application de la ther¬ 
modynamique en temps fini afin d’optimiser les 
conditions du fonctionnement du moteur de Stir¬ 
ling endo- et exo-irreversible, en vue d’obtenir un 
regime de puissance maximale. 

Par une etude de sensibility, nous avons mis en 
evidence l’influence des pertes de chaleur et des 
irreversibilitys internes dans le cas endoryversible, 
sur les indices de performance du moteur, c’est-a- 
dire la puissance mycanique foumie P et sur le 
rendement thermique 77. Nous avons aussi suivi 
les effets de ces irreversibihtys sur_ la valeur 
maximale de la puissance mecanique (P ma x) et sur 
le rendement correspondant r] opt . 

Les resultats obtenus sur la base du modele 
propose mettent en evidence le fait que, dans le 
cas endo-irreversible, pour le regime de puissance, 
les indices de performance du moteur de Stirling 
baissent d’environ 40 % par rapport au cas en- 
doreversible. II resulte aussi que, des trois sources 
d’irreversibilitys internes consideryes (cycle non 
quasi statique, conductance interne de la machine 
et regenyrateur non adiabatique), la diminution 
la plus importante des performances (« 30 %) est 
due au rygynerateur non adiabatique. Dans le cas 
du regime de puissance maximale, les resultats du 
modele endo- et exo-irreversible mettent en evidence 
les sources de non-concordance entre les resultats 
theoriques et les resultats experimentaux obtenus 
pour y opt , dans le cas du modele endoreversible 

[10]. Tous ces resultats sont en concordance avec 
les donnees experimentales connues pour le moteur 
de Stirling [1-4]. 
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